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THE  INTERACTION  OF  STATISTICS  AND  GEOLOGY  --  FINITE  DEFORMATIONS 
by  G.S.  Watson,  Princeton  University 

ABSTRACT 

The  paper  describes  briefly  the  Instances  where  the  analysis 
of  geological  data  has  required  the  development  of  new  statistical 
theory  and  methods.  A  new  instance  of  this  fruitful  interaction 
of  statistics  and  geology  is  then  given.  §2  develops  the  theory 
of  finite  deformations .  Some  of  the  classical  and  novel  statistical 
problems  which  arise  when  objects  embedded  in  a  deformed  rock  are 
measured  for  strain  estimation  are  discussed  in  §3. 
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1.  INTRODUCTION 

Sir  Ronald  Fisher  (1953a)  in  a  Presidential  address  to  the 
Royal  Statistical  Society  on  the  "Expansion  of  Statistics"  used 
as  his  first  example  the  statistical  thinking  of  Lyell  whose 
"Principles  of  Geology"  is  now  150  years  old  and  whose  publication 
is  celebrated  In  this  volume.  While  this  is  an  example  of  statis¬ 
tics  helping  Geology,  it  did  not  lead  to  developments  in  Statistics 
In  this  paper  we  wish  to  consider  (with  key  references)  some  In¬ 
stances  of  mutually  beneficial  Interactions  and  to  describe  some 
specific  results  of  Interest  to  structural  geologists  and 
statisticians. 

In  the  same  year  Fisher  published  his  famous  paper  on  "Disper¬ 
sion  on  a  Sphere"  (1953b).  Its  primary  purpose  was  to  provide 
statistical  methods  for  the  new  data  of  palaeomagnetism  but  he 
also  used  It  to  Illustrate  his  theory  of  Fiducial  Inference. 

Though  orientation  data  are  not  peculiar  to  Geology  and  Geophysics, 
it  was  these  subjects  which  first  demanded  such  statistical  methods 
This  led  to  a  new  area  of  statistical  theory  and  practice  and  to 
the  improvement  of  data  analysis  in  Geology.  This  story  is  now 
quite  well  known  --  see  e.g.  the  paper  of  Watson  (1970)  and  the 
book  by  Mardia  (1972).  Watson's  paper  contains  brief  sections  on 
orientation  analysis  in  Palaeomagnetism,  Sedlmentology ,  Structural 
Geology  and  Petrofabrics .  An  update  of  this  would  discuss  among 
other  things  the  relation  to  Plate  Tectonics  and  the  use  of  orien¬ 
tation  statistics  to  estimate  the  angular  momentum  vectors  of 


the  plates  and  their  finite  displacements,  (see  e.g.  Molnar  et  al  ( 

The  mutually  beneficial  Interaction  of  Statistics  and  Geology 
and  Mining  In  the  description  of  the  fabric  of  sediments  and  the 
dispersion  of  commercially  valuable  ores  Is  also  well  known  and 
particularly  associated  with  Matheron  and  the  Fontainebleau  School. 

Less  well  known  perhaps  Is  the  enormous  Interaction  In  quan¬ 
titative  Seismology,  mainly  due  to  the  towering  figure  of  Sir 
Harold  Jeffreys.  In  his  book  "The  Earth"  (1962-first  published 
1924)  he  remarks  (p.395):  "It  Is  astonishing  that  experimenters 
will  spend  months  In  making  a  series  of  observations  and  grudge 
the  day  or  so  needed  to  present  the  results  In  an  Intelligible 
form".  This  led  him  to  write  two  books  "Scientific  Inference" 
(1931)  and  "The  Theory  of  Probability"  (1939).  The  latter  was  the 
first  large  scale  account  of  Statistics  from  a  Bayesian  point  of 
view.  Jeffrey's  did  not,  as  a  classical  scientist,  espouse  sub¬ 
jective  priors.  To  obtain  his  priors  he  required  that  they  be  In¬ 
variant  under  functional  transformations  of  the  parameter  In 
question.  This  rule  Is  the  basis  of  some  modern  accounts  of 
Bayesian  methods  e.g.  Box  and  Tlao  (1973). 

But  there  Is  much  more  In  the  1939  book  that  Is  of  current 
Interest  to  statisticians  whatever  their  views  on  Inference.  For 
example,  he  remarks  that  while  different  observers  looking  at  a 
seismogram  will  usually  agree  to  a  second  or  so  when  a  particular 
wave  train  arrived,  occasionally  there  are  discrepancies  of  10  or 
so  seconds.  Thus  the  normal  distribution  for  arrival  times  will 
not  be  followed.  He  models  this  both  by  a  distribution  with 


heavier  tails  than  normal  and  as  a  "contaminated"  distribution. 

He  points  out  that  the  maximum  likelihood  (m.l.)  estimator  G 
for  e  In  the  density  f(x-e)  satisfies,  with  data  Xj,...,xn, 
Ef'(x.|-8)/f(Xj-$)  ■  0  .  He  writes 

»  w(x-e) 

so  that  the  m.l.  equation  can  be  written  as 
E(x-§)w(x-e)  -  0  , 

and  that  0  Is  thus  a  weighted  mean  of  the  observations.  He 
goes  on  to  regard  the  observations  as  a  mixture  of  Gaussian  and 
long-tailed  distributions.  These  Ideas  are  precisely  those  sug¬ 
gested  much  later  by  Tukey,  Huber  and  others  concerned  with  Robust¬ 
ness  --  see  e.g.  Andrews  et  al.  (1972),  Hosteller  and  Tukey  (1977). 

Modern  theory  and  data  analytical  methods  In  Quantitative 
Seismology  may  be  best  studied  now  In  the  two-volume  work  of  Ak 1 
and  Richards  (1980).  In  an  early  paper  however,  Scheldegger  (1964) 
suggested  the  method  for  estimating  the  fault  plane  movement  along 
which  caused  the  earthquake  being  analyzed.  In  this  paper  he  In¬ 
troduced  the  probability  density  proportional  to  exp<cos2e  which 
was  Independently  suggested  by  Dlmroth  (1963)  and  Watson  (1965) 
for  other  orientation  problems.  Scheldegger  wrote  other  papers 
on  this  area  with  Fara  with  whom  he  also  wrote  the  first  paper  on 
the  stochastic  description  of  sedimentary  fabr1cs--Fara  and  Scheldegger 
(1961).  However  most  selsmologlcal  data  analysis  Is  either 
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generallzed  regression  analysis  or,  more  Interestingly,  very 
elaborate  time  series  analysis.  References  to  the  latter  topic 
may  also  be  found  Brllllnger  (1975)  --  the  review  of  the  uses  of 
spectral  analysis  In  Geophysics  by  Tukey  (1965)  should  be 
especially  noted.  Seismic  methods  for  exploration  raise  many 
problems  --  see  e.g.  the  papers  and  books  by  Enders  Robinson  e.g. 
(1978). 

The  analysis  of  seismic  signals  may  be  pursued  to  Investigate 
their  source  (e.g.  earthquake  or  explosion)  or  to  Investigate  the 
media  they  pass  through  from  source  to  receivers.  This  theory  Is 
of  course  based  on  the  theory  of  elastic  deformations.  Structural 
geologists  also  use  the  theory  of  elasticity  to  unravel  the  history 
of  rock  formations. 

In  the  next  section  a  brief  account  Is  given  of  a  simple  form 
of  elastic  deformation,  homogeneous  strain.  In  section  three, 
some  statistical  problems  suggested  by  homogeneous  strain  are 
described  as  another  example  of  a  fruitful  Interaction  between 
Statistics  and  Geology. 


2.  HOMOGENEOUS  STRAIN 


While  In  seismology  one  usually  deals  with  small  reversible 
deformations,  over  long  time  periods  rocks  acquire  large  strains 
--  see  e.g.  Ramsey  (1967).  These  might  be  the  result  of  a 
sequence  of  linear  deformations  or  homogeneous  strains.  In  this 
section  we  summarize  the  description  of  such  strains  In  p 
dimensions.  In  practice  p*2  or  3  . 

A  physical  deformation.  If  linear,  Is  characterized  by 
matrix  £  with  whose  determinant  |£|  >  0  .  For  the  transforma¬ 
tion  must  grow  from  %  and  be  1  -  1  .  Since  Its  determinant 
cannot  change  sign  and  Is  Initially  positive.  It  must  remain  so. 
Any  two  points  Initially  joined  by  a  vector  of  length  sQ 

are  separated  by  s  after  the  deformation,  where 

■  Vio 

— 2 —  — n -  d) 

When  the  deformation  Is  small,  the  left  hand  side  of  (l)  Is 
approximately  2(s-Sq)/sq  .  Hence  the  strain  tensor  or  matrix 
j[  Is  defined  by 

£  .  l/2(£'£  -  l)  (2) 

If  the  eigen  vectors  and  values  of  £  are  and 
respectively  (1  ■  l...p)  ,  the  are  called  the  extensions 
since  they  equal  (s-Sq)sq  when  £q  Is  In  the  direction  . 
The  n1  may  be  negative  If  there  Is  compression  In  the  direc¬ 
tion  ,  positive  If  there  Is  tension  and  >  -1/2  . 
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In  spaces  of  odd  dimension,  since  |£f  >  0,  £  must  have 
at  least  one  positive  eigenvalue  and  to  each  positive  eigenvalue 
corresponds  an  Invariant  direction.  There  are,  however,  p 
mutually  orthogonal  directions  which  remain  so  after  deformation 
and  It  Is  easy  to  verify  that  these  are  defined  by  the  eigen¬ 
vectors  of  £  ,  the  principal  strain  axes.  Deformation 

corresponds  to  a  rotation  of  these  axes  and  extensions  along 
them. 

A  strain  £  Is  said  to  be  pure  If  and  only  if  It  leaves 
Invariant  p  orthogonal  directions  say.  Thus  these 

vectors  must  be  right  eigenvectors  of  £  with  non-zero 

roots,  ^.....cy  say.  Writing  JJ  *  [j^ . j£p],  JJ(a)  • 

d1ag(aj , . . . ,ap) ,  we  have  the  symmetric  form  for  £ 

H  ■  )«(•»)«'  •  <3> 

In  this  case, 

l  -  l/2<W<a2>U'  -  Ip)  <4> 

The  general  deformation  |  may  now  be  decomposed  Into  a 
rotation  *  1  »  followed  by  a  pure  strain  ^  ,  or  a 

pure  strain  £2  followed  by  a  rotation  J^*  |$2I  "  *  •  These 
are  just  the  polar  decompositions  of  £  .  Since  £  ■  «  £2£2  » 

*  $2  "  +  Ip  so  that  8  suitable  choice  for  £2  Is 

(5) 

j 


(6) 


where  s1  *  /2 n i  + 1  .  Then  Is  given  by 

h  ■ 

The  other  factorization  uses  the  eigenvectors  of  ^ 

Deformations  are  often  described  in  terms  of  their  effects 
on  a  quadric  surface,  -  c2  centered  at  £  .  It  will  be 

moved  and  become  the  surface 

t'llU'Vl'V  * c2  • 

Thus  the  Image  of  a  sphere,  called  the  strain  ellipsoid,  is 

*  c2  ,  (7) 

and  the  ellipsoid  that  has  a  spherical  image,  called  the 
reciprocal  strain  ellipsoid.  Is 

*  c"  •  (8) 

The  principal  axes  of  (7)  and  (8)  are  those  of  ^  and  £2 
respectively. 

If  the  strain  Is  small,  £  •  |p  ♦  £,  £  small,  so  that, 
approximately, 

•  (9) 

It  Is  customary,  when  dealing  with  small  strains,  to  write  the 
displacement  ^  *  8x  •  x  of  a  point  x  as 
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~~l  £  • 


*  &  +  u  » 


do) 


where  the  skew-symmetric  matrix  $  describes  the  rotation  and 
£  the  pure  strain,  which  are  now  additive.  Thus,  approximately, 

$1  *  $2  *  Hi  +  -Ip  •  $1  *  $2  *  £  +  -Ip 


(ID 
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3.  STATISTICAL  PROBLEMS  WITH  HOMOGENEOUS  STRAINS 

Strain  must  usually  in  Geology  be  determined  by  the  deforma¬ 
tion  of  objects  embedded  in  the  rock.  It  is  e.g.  rarely  possible 
to  use  devices  like  strain  gauges.  Further  the  rock  is  opaque 
so  we  are  usually  restricted  to  data  obtainable  from  the  surface 
of  plane  sections  cut  through  rock  samples.  We  must  of  course 
know  something  of  the  original  shape  and  distribution  of  the 
objects . 

Point-like  objects  distributed  at  random  (by  a  Poisson 
process)  with  unknown  mean  density  X  per  unit  volume  in  the 
rock  would  be  useless.  For  the  points  in  any  region  go  with  it 
under  the  transformation  £  while  the  volume  of  the  region  in¬ 
creases  by  the  factor  |J£|  so  the  transformed  field  of  points 
is  again  Poisson  but  with  density  X/|£|  .  If  X  were  known  |£| 
could  of  course  be  found.  |J£|  is  called  the  dilatation  of  the 
strain.  To  get  more  information  about  £  ,  the  points  distribution 
must  have  more  structure,  as  will  be  seen  below.  But  interesting 
mathematical  problems  remain  In  the  study  of  affine  transformations 
of  stationary  point  processes,  more  general  than  the  Poisson 
process . 

If  we  could  measure  the  initial  positions  ^  at  n  points 
and  then  measure  the  final  positions  ,  jji  may  be  estimated 
when  e.g.  the  measurement  errors  In  both  cases  are  Independently 
Gaussian  with  mean  vectors  zero  and  covariance  matrices  o  I  -- 
see  Gleser  and  Watson  (1973).  Unfortunately  though  this  might 
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be  used  in  glaciology,  it  is  not  of  use  with  ancient  rocks.  This 

problem  is  a  generalization  of  the  "errors  in  variables"  problem 

in  a  situation  where  it  admits  a  clean  solution.  The  simple 

form  (p=l)  of  this  problem  may  be  found  in  textbooks  (see  e.g. 

Theil  (1971))  on  Econometrics:  y=B£+f.  x=£+e  where  the  errors 

of  measurement  e  and  f  of  x  and  y  are  independent  Gaus- 

2  2  2 

sians  with  means  zero  and  variances  and  o  .  If  a  and 

x  y  x 

2 

ay  and  £  are  arbitrary  and  unknown,  there  is  no  satisfactory 

solution  basically  because  each  new  observation  introduces  a  new 

unknown  £  .  But  if  a *o  ,  there  is  a  satisfactory  solution. 

x  y 

Useful  information  about  the  strain  matrix  may  however  be 
obtained  f rom poi nts  wi th  unknown  initial  positions  but  arranged 
in  special  patterns  e.g. 

(i)  lines  of  known  initial  length 
(11)  angles  of  known  size 
(ill)  spherical  or  nearly  spherical  bodies. 

As  an  example  of  (i)  ,  consider  a  crystal  of  tourmaline  embedded 
in  a  rock  which  is  then  stretched  in  the  direction  of  this  long 
narrow  crystal.  The  crystal  will  break  into  a  number  of  pieces 
If  the  strain  is  large  enough.  The  initial  length  is  the  sum  of 
the  pieces.  The  final  length  Is  the  distance  from  the  beginning 
of  the  first  fragment  to  the  end  of  the  last.  An  example  of  (ii) 
would  be  a  brachiopod  whose  shape  Is  roughly  like  a  semi-circle 
before  deformation.  The  radius  perpendicular  to  the  base  Is, 
after  deformation,  inclined  at  a  measurable  angle  to  the  base. 


A 
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lt  is  estimated  that  70%  of  all  strain  studies  are  of  type  (iii). 
Oolites  (concretory  structures)  are  a  common  source  of  data. 

They  are  initially  roughly  spherical  so  their  final  shape  will 
be  (roughly)  the  strain  eillipsoid.  Cutting  plane  sections 
through  solids  and  studying  the  statistics  of  the  intersections 
with  the  embedded  bodies  leads  to  well  known  problems.  The 
theory  of  such  methods  usually  assumes  purely  random  planes-- 
however  structural  geologists  usually  make  cuts  related  to  the 
roughly  known  strain  axes  which  changes  the  methods. 

Statistical  methods  for  these  problems  were  outlined  in 
Watson  (1968)  and  will  appear  in  Watson  (1981).  Many  papers 
have  appeared  in  the  interim  --  see  e.g.  the  discussion  organized 
by  Ramsey  and  Wood  (1976)  and  the  references  therein.  One  paper 
--  Owens  (1973)  --  is  particularly  interesting  and  a  sequel  to 
March  (1932).  It  discusses  the  modification  under  strain  of  an 
angular  density  distribution. 

Let  f0(£)  denote  the  density  of  lines  with  direction  £  , 

1  in  the  solid  before  deformation.  Then  f0(4)<$w0 
proportional  to  the  number  of  lines  In  the  cone  about  £  with 
solid  angle  6«0  .  After  deformation,  the  axis  of  this  cone  will 
be  parallel  to  and  Its  volume  will  be  |$|  times  its 

previous  value  6u0/3  .  If  the  new  cone  has  solid  angle  6u>t, 

Its  length  Is  ||^||  ,  and  so  Its  volume  Is  || J^||  *6u)i/3  .  Hence 

ISUa,,  •  W«“. 

~r  t 
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o  r 


I«4I* 

ur 


6  co  i 


(12) 


The  new  density  fj  of  line  directions  is  then  determined  by 
the  fact  that  the  number  of  lines  in  each  cone  are  equal  i.e., 
by 

VM'IMD4"*  ■  f0U>6“o  (n) 

Using  (2.12). 

f !<«/■« I  '  VW’  (14) 

Iff- 

as  obtained  by  Owens  in  different  notation.  Hence  putting 

IE  s  WW  •  1  *  JfV*Bro*  • 

(14)  may  be  rewritten  as 

f i <n>)  ■  vrVifci  '  isiirVi  <«> 

A  full  discussion  of  this  method  of  modifying  spherical 
distributions  will  be  given  in  Watson  (1981).  Note  however  that 

if  f0(4)-f0(-4)  »  the  same  1S  true  ^1  and  that  ^  8  *s  a 

rotation  (14)  reduces  to  fj(B£)  *  fQ(£)  •  These  are  both 
necessary  checks.  Furthermore  if  fQ  *  (4ir)"*  ,  the  uniform 
distribution,  in  (15),  the  non-obvious  formula 

/nr1*  I  "*dw  -  |  B I  (16) 

m  •  i 

is  a  consequence  of  the  integral  of  fj(m)  being  unity. 


-13- 


Much  of  the  profit  that  statisticians  may  derive  from  the 
study  of  homogeneous  strain  comes  from  the  physical  intuition 
it  associates  with  linear  transformations.  To  those  interested 
in  orientation  or  directional  statistics  it  may  seem  awkward 
because  unit  vectors  before  strain  usually  do  not  have  unit  length 
afterwards.  However  many  directional  distributions  are  often  best 
thought  of  as  distributions  of  vectors  ^  ,  conditional  upon 
||  v  1 1  =  1  or  marginally  i.e.  after  integrating  out  their  length  and 
retaining  only  their  direction.  For  example,  the  Arnold-Fisher- 
von  Mises  distribution  with  density  proportional  to  expic^'^ 
where  k_>o,  ||£||  *  ||^||  -  1  can  be  obtained  exactly  as  the 
marginal  distribution  of  «  R£  ,  a  Gaussian  vector  with  mean  p  , 
covariance  matrix  a2£p  »  when  £  is  fixed.  But  it  may  be 
obtained  to  a  very  good  approximation  as  the  marginal  distribution 
of  l  --  the  distribution  actually  obtained  is  the  angular  Gaussian. 
For  mathematical  and  numerical  comparisons  see  Watson  (1980a, b). 

As  a  matter  of  history  it  is  interesting  to  note  that  Arnold's 
Thesis  (1941)  was  motivated  by  the  dissatisfaction  of  J.F.  Bell 
with  the  then  available  methods  for  analyzing  data  on  the  preferred 
direction  of  the  optical  axes  of  crystals  In  rocks.  Prof.  Bell 
subsequently  had  a  distinguished  career  at  Johns  Hopkins  in 
Mechanics.  Geology  has  certainly  benefltted  from  the  development 
of  these  methods.  Thus  this  Is  an  example  of  the  fruitful  Inter¬ 
action  of  Statistics  and  Geology. 


-14- 


REFERENCES 

AKI,  K.  ,  RICHARDS,  P.G.  (1980).  "Quantl tati ve  Seismology,  Theory 
and  Methods”,  Vols.  1,2,  W.H.  Freeman  &  Co.,  San  Francisco, 

Calif.  ,  pp.  932. 

ANDREWS,  D.F.,  BICKEl,  P.J.,  HAMPEL,  F.R.,  HUBER,  P.J.,  ROGERS,  W.H., 
TUKEY ,  J.W.  (1972).  "Robust  Estimates  of  Location",  Princeton 
University  Press,  Princeton,  N.J.,  pp.  369. 

ARNOLD,  K.J.  (1941).  "On  Spherical  Probability  Distributions", 

Ph.D.  Thesis,  Mass.  Inst.  Technology.,  pp.  42. 

BOX,  G.  and  TIAO,  G.C.  (1973).  "Bayesian  Inference  In  Statistical 
Analysis",  Addl son-Wesl ey ,  Reading,  Mass.,  pp.  588. 

BRILLINGER,  D.R.  (1975).  "Time  Series:  Data  Analysis  and  Theory", 

Holt,  Rinehart  &  Winston,  Inc.,  New  York.,  pp.  500. 

DIMROTH,  E.  (1963).  " Fortschl tte  de  Gefiigestati  stl  k" ,  N.  Jb.  Mlnpr. 

Mh. ,  Heft  8,  pp.  186-192. 

FARA,  H.D.,  SCHEIDEGGER,  A. E.  (1961).  "Statisical  Geometry  of 

Porous  Media",  Jour.  Geoohys.  Res..  Vol.  66,  No.  10,  pp.  3279-3284. 

FISHER,  SIR  RONALD  A.  (1953a).  "Dispersion  on  a  Sphere",  Proc. 

Roy.  Soc.  Lond..  Vol.  A217,  pp.  295-305. 

FISHER,  SIR  RONALD  A.  (1953b).  "The  Expansion  of  Statistics", 

J.  Roy.  Stat.  Soc..  Vol.  116,  Pt.  I,  pp.  1-9. 

GLESER,  L.J.,  WATSON,  G.S.  (1973).  "The  Estimation  of  a  Linear 
Transformation",  Blometrlka ,  Vol.  60,  Pt.  3,  pp.  524-535. 

JEFFREYS,  SIR  HAROLD  (1931).  "Cartesian  Tensors",  Cambridge 
University  Press.,  pp.  92. 


-15- 


JEFFREYS,  SIR  HAROLD  (1937).  "Scientific  Inference",  Cambridge 
University  Press.,  pp.  103. 

JEFFREYS,  SIR  HAROLD  (1939).  "Theory  of  Probability",  Clardendon 
Press,  Oxford.,  pp.  380. 

MARCH,  A.  (1932).  "Mathemati sche  Theorle  der  Regelung  nach  der 

Korngestalt  be  Afflnen  deformation",  Z.  Krlst.,  Vol .  81,  Heft  3-4, 
pp. 285-297. 

MARDIA,  K.V.  (1972).  "Statistics  of  Directional  Data",  Academic 
Press,  London  and  New  York.,  pp.  357. 

MINSTER,  J.B.,  JORDAN,  T.H.,  MOLNAR,  P.,  HAINES,  E.  (1974). 

"Numerical  Modelling  of  Instantaneous  Plate  Tectonics", 

Geophysics  J.R.  Astr.  Soc.  (1974),  pp.  541-76.  No.  3. 

MOSTELLER,  F.,  TUKEY,  J.W.  (1977).  "Data  Analysis  and  Regression", 
Addlson-Wesley,  Reading,  Mass.,  pp.  588. 

NYE,  J.F.  (1957).  "Physical  Properties  of  Crystals",  Clardendon 
Press,  Oxford.,  pp.  322. 

OWENS,  W.H.  (1973).  "Strain  Modification  of  Angular  Density 

Distributions",  Tectonphyslcs ,  Vol.  16,  No.  3-4,  pp.  249-261. 

RAMSEY,  J.G.  (1967).  "The  Folding  and  Fracturing  of  Rocks", 

McGraw-Hill  Book  Co.,  New  York.,  pp.  568. 

RAMSEY,  J.G.  and  WOOD,  D.S.  (1976).  "A  Discussion  of  Natural 

Strain  and  Geological  Structure",  Phil.  Trans.  Roy.  Soc.  Lond.. 
Vol.  283,  pp. 1-344. 

ROBINSON,  E.A.,  SILVIA,  M.T.  (1978).  "Digital  Signal  Processing 
and  Time  Series  Analysis",  Holden-Day,  Inc.,  San  Francisco, 

Calif.,  pp.  405. 


SCHEIDEGGER,  A.E.  (1964).  "The  tectonic  stress  and  tectonic 
motion  direction  In  Europe  and  Western  Asia  as  calculated 
from  earthquake  fault  plane  motions".  Bull.  Sels.  Soc.  Amer. , 

Vol .  54,  No.  5,  Part  A,  pp. 1519-1528. 

THEIL ,  H.  (1971).  "Principles  of  Econometrics",  J.  Wiley  &  Sons, 

New  York.  ,  pp.  736. 

TUKEY ,  J.W.  (1965).  "The  Uses  of  Numerical  Spectrum  Analysis  In 
Geophysics",  Bull.  1 . S . I .  ,V  35e .  Session,  pp. 267-307. 

WATSON,  G.S.  (1965).  "Equatorial  Distributions  on  a  Sphere", 

Blometrl ka ,  Vol.  52,  Pts.  142,  pp.  193-201. 

WATSON,  G.S.  (1968).  "Orientation  Statistics  In  the  Earth  Sciences", 
Technical  Report  No.  93,  The  Department  of  Statistics, 

The  Johns  Hopkins  University,  97  pp  6  6  figures. 

WATSON,  G.S.  (1970).  "Orientation  Statistics  In  the  Earth  Sciences", 

Bull,  Geol .  Inst.,  University  of  Uppsala, New  Series  2:9,  pp.  73-89. 

WATSON,  G.S.  (1980a).  "The  Estimation  of  Palaeomagnetlc  Pole 
Positions",  Technical  Report  No.  169,  Series  2,  July  1980 
(To  appear  In  volume  In  honour  of  C.R.  Rao's  60th  birthday  14  pp.) 

WATSON,  G.S.  (1980b).  "Distributions  on  the  Circle  and  Sphere", 
Technical  Report  No.  176,  Series  2,  November  1980  (To  appear 
In  volume  In  honour  of  P.A.P.  Moran's  65th  birthday  37  pp.) 

WATSON,  G.S.  (1981).  "Orientation  Problems" 

(To  be  read  to  Royal  Statistical  Society  In  October  1981  and 
to  appear  J.  Royal  Statistical  Society  Series  B). 


UNCLASSIFIED 


AECUftITv  CLASSIFICATION  OF  THIS  RAGE  (Whu  Data  EMafaG) 


REPORT  DOCUMENTATION  PAGE 


A.  TlTLC  (m*  tuMlUJ 


THE  INTERACTION  OF  STATISTICS  AND  GEOLOGY- 
FINITE  DEFORMATIONS 


AUTNORfa) 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


s.  reorient**  catalog  number 


Technical 


Geoffrey  S.  Watson 


N00014-79-C-0322 


Department  of  Statistics  * 
Princeton  University 
Princeton,  NJ  08544 


n.  CONTROLLING  OFFlCC  NAME  ANO  ADDRESS 

Office  of  Naval  Research  (Code  436) 
Arlington,  VA  22217 


MONITORING  AGENCY  NAME  A  ADORES*’!!  GJ /tonal  tnm  CaMrafURf  OMm) 


IS.  RERORT  DATE 

November  1980 


ts.  NUMBER  OF  RAGES 

Sixteen! 16) 

CURITy  CLASS.  fa# 1 

Unclassified 


APPROVED  FOR  PUBLIC  RELEASE:  DISTRIBUTION  UNLIMITED. 


17.  DISTRIBUTION  STATEMENT  fa#  *»  aAairacI  anfaraG  <a  Mask  SB,  If  mUanrnt  I 


IS-  SURRLSMENTARV  NOTES 


To  appear  in  a  volume  celebrating  the  150th  anniversary  of  the  publication 
of  Lyell’s  book  on  "Stratigraphy"  which  Inspired  Charles  Darwin. 


ft.  KEY  WORM  f Centime*  m  revere*  eMi  If  neeeeeerf  me*  JAmtfff  if  llteft  MBler) 

palaeomagnetism,  Harold  Jeffreys  text,  plate  techtonics,  homogeneous  strain, 
structural  geology. 


ABSTRACT  (Cmttmm  m  revere*  ef#t  tfneteeif 


The  paper  describes  briefly  the  Instances  where  the  analysis  of  geological 
data  has  required  the  development  of  new  statistical  theory  and  methods. 

A  new  Instance  of  this  fruitful  Interaction  of  statistics  jnd  geology  Is 
then  given v— |S  develops  |he  theory  of  finite  deformations.  Some  of  the 
classical  and  novel  statistical  problems  which  arise  when  ‘objects  embedded 
In  a  deformed  rock  are  measured  for  strain  estimation  are  discussed 


hi 


I  JAN  71 


Edition  of  inoveeie 
Ml  Sltt-LP<Sl44ESt 


UNCLASSIFIED 


